Wide binaries as a critical test for Gravity theories 
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Abstract. Assuming Newton's gravity and GR to be valid at all scales leads to the dark matter 
hypothesis as a requirement demanded by the observed dynamics and measured baryonic content 
at galactic and extragalactic scales. Alternatively, modified gravity scenarios where a change 
of regime appears at acceleration scales a < ao have been proposed. This modified regime at 
a < ao will generically be characterised by equilibrium velocities which become independent 
of distance. Here we identify a critical test in this debate and we propose its application to 
samples of wide binary stars. Since for IMq systems the acceleration drops below ao at scales 
of around 7000 AU, a statistical survey of wide binaries with relative velocities and separations 
reaching 10 4 AU and beyond should prove useful to the above debate. We apply the proposed 
test to the best currently available data. Results show a constant upper limit to the relative 
velocities in wide binaries which is independent of separation for over three orders of magnitude, 
in analogy with galactic flat rotation curves in the same a < ao acceleration regime. Our results 
are suggestive of a breakdown of Kepler's third law beyond a « ao scales, in accordance with 
generic predictions of modified gravity theories designed not to require any dark matter at 
galactic scales and beyond. 



1. Introduction 

Over the past years the dominant explanation for the large mass to light ratios inferred for 
galactic and meta-galactic systems, that these are embedded within massive dark matter halos, 
has begun to be challenged. Direct detection of the dark matter particles, in spite of decades of 
extensive and dedicated searches, remains lacking. Numerous alternative theories of gravity have 
appeared (e.g. TeVeS of [1], F(R) theories e.g. [8]) mostly grounded on geometrical extensions 
to General Relativity, and leading to laws of gravity which in the low acceleration regime, 
mimic the MOdified Newtonian Dynamics (MOND) fitting formulas e.g. [7J. Similarly, [6] have 
explored MOND not as a modification to Newton's second law, but as a modified gravitational 
force law in the Newtonian regime, finding a good agreement with observed dynamics across 
galactic scales without requiring dark matter. Recently [2] have constructed an f(R) extension 
to general relativity which in the low velocity limit converges to the above approach. 

A generic feature of all such modified gravity schemes is the appearance of an acceleration 
scale, ao, above which classical gravity is recovered, and below which the dark matter mimicking 
regime appears. This results in a general prediction; all systems where a >> ao should appear as 
devoid of dark matter, and all systems where a « ao should appear as dark matter dominated. 
It is interesting that no a >> ao system has ever been detected where dark matter needs to 
be invoked. On the other hand, the latter condition furnishes a testable prediction in relation 
to the orbits of wide binaries. For test particles in orbit around a IMq star, in circular orbits 



of radius s, the acceleration is expected to drop below ao ~ 1.2 x W~ 10 m/s 2 for s >7000 
AU= 3.4 x 10 _2 pc. This provides a test for the dark matter/modified theories of gravity debate; 
the relative velocities of components of binary stars with large physical separations should 
deviate from Kepler's third law. Seen as an equivalent Newtonian force law, beyond s ~7000 
AU the gravitational force should gradually switch from the classical form of Fjy = GM/s 2 
to Fmg = (GMao) 1 / 2 /s, and the orbital velocity, V 2 /s = F, should no longer decrease with 
separation, but settle at a constant value, dependent only on the total mass of the system through 
V = (GMao) 1 / 4 . Under modified gravity theories, binary stars with separations beyond 7000 
AU should exhibit "flat rotation curves" and a "Tully-Fisher relation", as galactic systems in 
the same acceleration regime do. 

We apply this test to the binaries of two recent catalogues containing relative velocities and 
separations. The two catalogues are entirely independent in their approaches. The first one, jS], 
comprises 280 systems from the Hipparcos satellite, while the second, [3], identifies 1,250 ones 
from the Sloan Digital Sky Survey (SDSS). Details can be found in 

2. Relative velocity distributions for wide binaries under Newtonian expectations 

The Newtonian prediction for the relative velocities of the two components of binaries having 
circular orbits, when plotted against the binary physical separation, s, is for a scaling of 
AU oc s~ 1//2 , essentially following Kepler's third law, provided the range of masses involved 
were narrow. In a relative proper motion sample however, only two components of the relative 
velocity appear, as velocity along the line of sight to the binary leads to no proper motion. 
Thus, orbital projection plays a part, with systems having orbital planes along the line of sight 
sometimes appearing as having no relative proper motions. A further effect comes from any 
degree of orbital ellipticity present; it is hence clear that the trend for AU oc s" 1 / 2 described 
above, will only provide an upper limit to the distribution of projected AU vs. s expected in 
any real observed sample, even if only a narrow range of masses is included. One should expect 
a range of measured values of projected AU at a fixed observed projected s, all extending below 



Further, over time, the orbital parameters of binaries will evolve due to the effects of Galactic 
tidal forces and dynamical encounters with other stars in the field, specially in the case of wide 
binaries. To first order, one would expect little evolution for binaries tighter than the tidal 
limit of 1.7pc, and the eventual dissolution of wider systems. A very detailed study of all these 
points has recently appeared, These authors numerically follow populations of 50,000 IMq 
binaries in the Galactic environment, accounting for the evolution of the orbital parameters of 
each due to the cumulative effects of the Galactic tidal field at the Solar radius. Also, the effects 
of close and long range encounters with other stars in the field are carefully included, to yield 
a present day distribution of separations and relative velocities for an extensive population of 
wide binaries, under Newtonian Gravity. 

It is found that when many wide binaries cross their Jacobi radius, the two components 
remain fairly close by in both coordinate and velocity space. Thus, in any real wide binary 
search a number of wide pairs with separations larger than their Jacobi radii will appear. Finally, 
[5] obtain the RMS one-dimensional relative velocity difference, AVid, projected along a line 
of sight, for the entire populations of binaries dynamically evolved over 10 Gyr to today, as 
plotted against the projected separation on the sky for each pair. The expected Keplerian fall 
of AUid oc s -1 / 2 for separations below 1.7 pc is obtained, followed by a slight rise in AVid as 
wide systems cross the Jacobi radius threshold. AUid then settles at RMS values of 0.1km/ s. 
This represents the best currently available estimate of how relative velocities should scale 
with projected separations for binary stars (both bound and in the process of dissolving in the 
Galactic tides) under Newtonian gravity. We see that all we need is a large sample of relative 



the Newtonian limit, which for equal mass binaries in circular orbits gives AUv = 2 




proper motion and binary separation measurements to test the Newtonian prediction for the 
RMS values of the 1 dimensional relative velocities of [5], and to contrast the AV/v oc s -1 / 2 and 
the AVmg = cte- predictions for the upper envelope of the AV vs. s distributions. 

3. Observed wide binary Samples 

In the [8] catalogue wide binaries are identified by assigning a probability above chance alignment 
to the systems by carefully comparing to the underlying background (and its variations) in a 
5 dimensional parameter space of proper motions and spatial positions. We keep only binaries 
with a probability of non-chance alignment greater than 0.9. The binary search criteria used by 
the authors requires that the proposed binary should have no near neighbours; the projected 
separation between the two components is thus always many times smaller than the typical 
interstellar separation. We use the reported distances to the primaries, where errors are smallest, 
to calculate projected AV and projected s from the measured Afi and A9 values reported. 
Although the use of Hipparcos measurements guarantees the best available quality in the data, 
we have also removed all binaries for which the final signal to noise ratio in the relative velocities 
was lower than 0.3. 

We are left with a sample of 280 binaries, having distances to the Sun within 6 < d < 100 
in pc. The data show a perfectly flat upper envelope in a AV vs. projected s, [I]. The average 
signal to noise ratio for the data is 1.7, with an average error on AV of 0.83 km/s, which 
considering a 2a factor from the top of the distribution to the real underlying upper limit for 
the sample, results in 3 km/s as our estimate of the actual physical upper limit in AV. 

The Sloan low mass wide pairs catalogue (SLoWPoKES) of [3 J contains a little over 1,200 
wide binaries with relative proper motions for each pair, distances and angular separations. 
Also, extreme care was taken to include only physical binaries, with a full galactic population 
model used to exclude chance alignment stars using galactic coordinates and galactic velocities, 
resulting in an estimate of fewer than 2% of false positives. This yields only isolated binaries with 
no neighbours within many times the internal binary separation. Again, we use the reported 
distances to the primaries to calculate projected AV and projected s from the measured A/x, 
AO and d values reported by [3], to obtain a sample of 417 binaries. 

The upper envelope of the distribution of AV from this catalogue does not comply with 
Kepler's third law. As was the case with the Hipparcos sample, the upper envelope describes a 
flat line, as expected under modified gravity schemes. The average signal to noise in AV for the 
[3J catalogue is 0.48, with an average error on AV of 12 km/s, which considering a 2a factor 
from the top of the complete distribution to the real underlying upper limit gives the same 
3km /s as obtained for the [8] Hipparcos catalogue. 

Figure (1) shows the RMS value of the one-dimensional relative velocity difference for both 
of the samples discussed. The error bars give the error propagation on Afi and d. We construct 
AVid by considering only one coordinate of the two available from the relative motion on the 
plane of the sky. Thus, each binary can furnish two AVid measurements, which statistically 
should not introduce any bias. Indeed, using only A/i; or only Afjn, or both for each binary, yields 
the same mean values for the points shown. The small solid error bars result from considering 
an enlarged sample where each binary contributes two AV\b measurements, while the larger 
dotted ones come from considering each binary only once, and do not change if we consider only 
Ayn or only A//&. The series of small log(s) interval data are for the Hipparcos catalogue of [8], 
while the two broader crosses show results for the [3J SDSS sample. 

The solid curve is the Newtonian prediction of the full Galactic evolutionary model of [5] for 
binaries, both bound and in the process of dissolving. Note that the results of this simulation 
deviate from Kepler's law for s larger than the Newtonian Jacobi radius of 1.7pc, whereas the 
discrepancy with the observed samples also occurs at much smaller separations. Even considering 
the large error bars, where each binary contributes only one AV\d value, we see eight points lying 



beyond la, making the probability of consistency between this prediction and the observations 
of less than (0.272) 8 =3 x 10" 5 . 

We obtain a constant RMS value for AV±d of 1 km/s, in qualitative agreement with 
expectations from modified gravity schemes. The vertical line marks a = ao; we see the 
data departing from the Newtonian prediction outwards of this line, and not before. The two 
independent catalogues, each using different sets of selection criteria, each perhaps subject to its 
own independent systematics, are consistent with the same result, a constant horizontal upper 
envelope for the distribution of relative velocities on the plane of the sky at an intrinsic value of 
3 km/s±l km/s, extending over 3 orders of magnitude in s, with a constant RMS AV\d value 
consistent with 1 fcm/s±0.5 km/s. This supports the interpretation of the effect detected as 
the generic prediction of modified gravity theories. 
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Figure 1. The solid curve gives 
the RMS values for one dimensional 
projected relative velocities as a 
function of projected separations, 
for the dynamical modelling of large 
populations of wide binaries evolv- 
ing in the Galactic environment, 
taken from [5]. The same quantity 
for the data from the catalogues 
analysed is given by the points with 
error bars; those with narrow log(s) 
intervals being from the Hipparcos 
sample of [8], and those two with 
wide log(s) intervals coming from 
the SDSS sample of [3] 



4. Conclusions 

We identify a critical test in the classical gravity/modified gravity debate, using the relative 
velocities of wide binaries with separations in excess of 7000 AU, as these occupy the a < ao 
regime. We present a first application of this test using the best currently available data. Results 
show constant relative RMS velocities for the binary stars in question. This is quantitatively 
inconsistent with detailed predictions of Newtonian dynamical models for binaries evolving in 
the local galactic environment. Our results are qualitatively in accordance with generic modified 
gravity models which explain galactic dynamics in the absence of dark matter. 
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